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We reanalyze deep inelastic scattering data of BCDMS Collaboration by including proper cuts of ranges with 
large systematic errors. We perform also the fits of high statistic deep inelastic scattering data of BCDMS, SLAC, 
NM and BFP Collaborations taking the data separately and in combined way and find good agreement between 
these analyses. We extract the values of both the QCD coupling constant a s (M z ) up to NLO level and of the 
power corrections to the structure function Fi- 



1. Introduction 

The deep inelastic scattering (DIS) leptons on 
hadrons is the basical process to study the values 
of the parton distribution functions which are uni- 
versal (after choosing of factorization and renor- 
malization schemes) and can be used in other 
processes. The accuracy of the present data for 
deep inelastic structure functions (SF) reached 
the level at which the Q 2 -dependence of logarith- 
mic QCD-motivated terms and power-like ones 
may be studied separately (for a review, see the 
recent paper PQ) and references therein). 

In the present letter we sketch the results of 
our analysis (see also [?]) at the next-to- 
leading order of perturbative QCD for the most 
known DIS SF F 2 (x, Q 2 ) 1 taking into account 
experimental data g]-[7| of SLAC, NM, BCDMS 
and BFP Collaborations. We stress the power- 
like effects, so-called twist-4 (i.e. ~ 1/Q 2 ) con- 
tributions. To our purposes we represent the 
SF F 2 (x,Q 2 ) as the contribution of the leading 
twist part F^ 00 (x,Q 2 ) described by perturba- 
tive QCD, when the target mass corrections are 
taken into account, and the nonperturbative part 
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where hi{x) is magnitude of twist-4 terms. 

1 Here Q 2 = —q 2 and x = Q 2 /(2pq) are standard DIS vari- 
ables, where q and p are photon and hadron momentums, 
respectively. 



Contrary to standard fits (see |3] and refer- 
ences therein) when the direct numerical cal- 
culations based on Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi equation [H] are used to evalu- 
ate structure functions, we use the exact solution 
of DGLAP equation for the Mellin moments of 
SF F$ D (x,Q 2 ) and the subsequent reproduc- 
tion of F% (x,Q 2 ) at every needed Q 2 -value 
with help of the Jacobi Polynomial expansion 
method |9I10| 2 (see similar analyses at the NLO 
level |1UI14| and at the next-next-to-leading order 
(NNLO) level and above [T5|h 

2. Results of fits 

We have studied in [3] the Q 2 -evolution of DIS 
structure function F 2 fitting all modern fixed tar- 
get experimental data at x > 10~ 2 and Q 2 > 1 
GeV 2 . From the fits we have obtained the value 
of the normalization a s (M^) of QCD coupling 
constant. 

First of all, we have reanalyzed the BCDMS 
data cutting the range with large systematic er- 
rors. As it is possible to see in the value of 
a s (M|) rises strongly when the cuts of systemat- 
ics were incorporated. In another side, the value 
of a s (Mg) does not dependent on the concrete 
type of the cut within modern statistical errors. 



2 We note here that there is similar method II 1 1 , based on 
Bernstein polynomials. The method has been used in the 
analyses at the NLO level in 1121 and at the NNLO level 
in 
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Figure 1. The values of a s (M§) and x 2 at dif- 
ferent Q 2 -values of data cuts in the fits based on 
nonsinglet evolution. The black (white) points 
show the analyses of data without (with) twist-4 
contributions. Only statistical errors are shown. 

To verify the range of applicability of per- 
turbative QCD, we analyzed firstly the SLAC, 
BCDMS, NM and BFP data without a contri- 
bution of twist-4 terms, i.e. when F2 = F^® 00 . 
We did several fits using the cut Q 2 > Q 2 cut and 
increased the value Q 2 ut step by step. We ob- 
served good agreement of the fits with the data 
when Q 2 cut > 10 -j- 15 GeV 2 (sec the Figs. 1 and 
2). Later we added the twist-4 corrections and 
fitted the data with the standard cut Q 2 > 1 
GeV 2 . We have found very good agreement with 
the data. Moreover, the predictions for a s (M§) 
in both above procedures are very similar (see the 
Figs. 1 and 2). 

The results for a s (M§) coincide for the both 
types of analyses: ones, based on nonsinglet evo- 
lution, and ones, based on combined singlet and 
nonsinglet evolution. They have the following 



Figure 2. The values of a s (M§) and \ 2 a t dif- 
ferent Q 2 -values of data cutes in the fits based 
on combine singlet and nonsinglet evolution. All 
other notes are as in Fig. 1. 

form: 

a s {M 2 z ) = 0.1177 ±0.0024 (total), (2) 

where the symbol "total" marks the total experi- 
mental error, which contains the sum of statistical 
error, systematic one and error of normalization 
in quadratures. 

We would like to note that the result J2J is 
in good agreement with the analysis |16| of com- 
bined HI and BCDMS data, which has been given 
by HI Collaboration very recently. The result for 
a s (Mz) is in good agreement also with the av- 
erage value for coupling constant, presented in 
famous Altarelli and Bethke reviews |17|. 
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